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Abstract 

The Doppler peaks (Sacharov peaks) in the angular power spectrum of 
the cosmic microwave background anisotropies, are mainly due to coherent 
oscillations in the baryon radiation plasma before recombination. Here we 
present a calculation of the Doppler peaks for perturbations induced by global 
textures and cold dark matter. We find that the height of the first Doppler 
peak is smaller than in standard cold dark matter models, and that its position 
is shifted to £ ~ 350. We believe that our analysis can be easily extended to 
other types of global topological defects and general global scalar fields. 


PACS numbers: 98.80-k 98.80.Hw 98.80C 


Presently there are two main classes of models to explain the origin of large scale 
structure formation. Initial perturbations can either be due to quantum fluctuations 
of a scalar held during an inhationary era|I|, or they may be seeded by topological 
defects formed during a symmetry breaking phase transition in the early universe p| . 
The CMB anisotropies are a powerful tool to discriminate among these models by 
purely linear analysis. Usually CMB anisotropies are parameterized in terms of Cfs, 
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defined as the coefficients in the expansion of the angular correlation function 
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For scale invariant spectra of perturbations £{£ + 1)C£ is constant on large angular 
scales, say £ ~ 50. Both inflation and topological defect models lead to approxi¬ 
mately scale invariant spectra on large scales. 

Large scale CMB anisotropies are mainly caused by inhomogeneities in the space- 
time geometry via the Sachs-Wolfe (SW) effect^]. On smaller angular scales (0.1° ~ 
6 ~ 2°) the dominant contribution comes from coherent oscillations in the baryon- 
radiation plasma prior to recombination. On even smaller scales the anisotropies are 
damped due to the finite thickness of the recombination shell, as well as by photon 
diffusion during recombination (Silk damping). 

Disregarding Silk damping, gauge invariant linear perturbation analysis leads 
toM 
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where <F and 'L are quantities describing the perturbations in the geometry and V 
denotes the peculiar velocity of the baryon fluid with respect to the overall Friedmann 
expansion. specifies the intrinsic density fluctuation in the radiation fluid. 

There are several gauge invariant variables which describe density fluctuations; they 
all differ on super-horizon scales but coincide inside the horizon. Below we use 
another quantity, D^, for the radiation density fluctuation. The variables and 
Dr are defined in Eq. (II.5.28b) and Eq. (II.5.27b) of ref. p[ respectively, for an 
arbitrary matter component a. Here r stands for the coupled baryon radiation 
fluidQ. Since the coherent oscillations giving rise to the Doppler peaks act only on 

sub-horizon scales, the choice of this variable is irrelevant for our calculation. 

^ Actually Dr = {5pr + Spb)/{pr + Pb) is not quite the variable for the temperature fluctuation, 
ST/T = {1/A)5pr/Pr- A short calculation shows that Dr is about 5% smaller than SprjPr- 
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$ , and in Eq. (J^) determine the anisotropies on large angular scales^, and 
have been calculated for both inflation and defect models i. 0. i a- Generically, 
a scale invariant spectrum is predicted and thus the SW calculations yield mainly 
a normalization for the different models. On the other hand the amplitude of the 
Doppler peak, which most probably will be measured in the near future, might be 
an important discriminating tool between them. In this Letter we present a com¬ 
putation for the Doppler contribution from global topological defects; in particular 
we perform our analysis for vrs-defects, textures [|l^], in a universe dominated by 
cold dark matter (CDM). We believe that our main conclusion remains valid for all 
global defects. 

The Doppler contribution to the CMB anisotropies is approximately given by0 

Doppler ^ 

~ (^rec; ^rec) ^i^^reci tree) ' ^7 (2) 

where x^ec = x — nto- In the previous formula n denotes a direction in the sky and t 
is conformal time, with to and tree the present and recombination times, respectively. 
Due to the inaccuracies mentioned in the previous footnotes, Eq. (|^) tends to under¬ 
estimate the amplitude of the hrst Doppler peak by up to 30%. On the other hand, 
we neglect Silk damping of perturbations, which leads to a slight overestimation. We 

thus are on the safe side, if we postulate that approximation leads to an error of 

^One might think that leads just to coherent oscillations of the baryon radiation fluid, but 
this is not the case. Note that, e.g., for adiabatic CDM models without source term one can derive 
(l/4)Dg’^^ = —(5/3)d> on super-horizon scales. Since for CDM perturbations, $ = —and dt' ~ 0, 
the usual SW result 5T/T = (l/3)d>(xr.ec, Gee) is recovered. Neglecting Dg^\ the result would be 
2'!' and therefore wrong by a factor of 6! 

^In principle this Doppler term has to be added to the SW contribution. But the SW contribu¬ 
tion decays on subhorizon scales (like At horizon scales, especially the last term in Eq. (|^, 

the integrated Sachs-Wolfe (ISW) effect, can be important. At i = 200 it contributes about 30% 
to the angular power spectrum for standard adiabatic CDM. Neglecting it, slightly shifts the first 
Doppler peak to smaller angular scales, £ « 220, as we have found by testing our code for the 
standard adiabatic CDM model. Since we will find here, that the first Doppler peak is lower than 
in this model, the ISW contribution might be higher. However, as we shall see below, the peak is 
at ^ = 365. Therefore we expect a suppression by (365/200)^ « 3, so that the ISW contribution 
to the peak is probably not much higher. We also neglect the contribution of the neutrino fluc¬ 
tuations. But since even the dark matter fluctuations yield only about 20% of the gravitational 
potential, we expect the neutrino fluctuations, which for standard models contribute about 20%, 
to be considerably smaller. 
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less than about 30% in the amplitude of the hrst Doppler peak and overestimates 
the value i of its position by less than 10%. The nice feature of Eq. (2) is that 
we will need only one simple scalar component of the defect stress-energy tensor to 
evaluate it. 

To determine Dr and V at tree, we consider a two-fluid system: baryons plus 
radiation, which prior to recombination are tightly coupled, and CDM. The evolution 
of the perturbation variables in a flat background, D = 1, is described by[|| 




fcT + k^D 

l-\-W 
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V’ + ^Ve = fcT 

D'^-3w^Dr = {1 + w)[3^^ - 3^'- kVr - ^ k-\l + ^)Vr] 

where subscripts r and c denote the baryon-radiation plasma and CDM, respectively; 
D, V are density and velocity perturbations; w = Pr/Pr, = p'r/Pr and p = Pr + Pc- 

The only place where the seeds enter this system is through the potentials T and $. 
These potentials can be split into a part coming from standard matter and radiation, 
and a part due to the seeds, T = + d's and <h = ‘h(c,r) + where and 

are determined by the energy momentum tensor of the seeds. In this way, the seed 
source terms will arise below[^]. 

From Eqs. (^) we derive two second order equations for Dr and De, namely 
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where F = k^{127iGa^) ^ + pr{l + w) + pc and S denotes a source term, which in 
general is given by S' = 47rGa^(p + 3p)®“'^. In our case, where the seed is described 
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by a global scalar field 0, we have S = 87rG(0')^. From numerical simulations one 
finds that the average of over a shell of radius k, can be modeled during the 
matter dominated era by[y] 


(|^'|">(fc,«) + alkt^ 
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with T] denoting the symmetry breaking scale of the phase transition leading to 
texture formation. The parameters in (|^) are A ~ 3.3, a ~ —0.7/(27r) and (3 ~ 
0.7/(27r)^. On super-horizon scales, where the source term is important, this fit is 
accurate to about 10%. As we argue later, analytical estimates support this finding. 
On small scales the accuracy reduces to a factor of 2. By using this fitQ in the 
calculation of and Dc from Eqs. (^), we effectively neglect the time evolution 
of phases of (0')^; the incoherent evolution of these phases may smear out subsequent 
Doppler peaks[|^, but will not affect substantially the height of the first peak. 

From Dr and D' we calculate the Doppler contribution to the O^’s according to 
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where ji denotes the spherical Bessel function of order i and stands for its deriva¬ 
tive with respect to the argument. The angular power spectrum + 1)C(^ yields 
the Doppler peaks. 

In order to solve Eqs. (^, (^) we need to specify initial conditions. For a given 
scale k we choose the initial time tin such that the perturbation is super-horizon and 
the universe is radiation dominated. In this limit the evolution equations reduce to 


D"-^D = 
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(9) 


'*Our fit is not valid in the radiation dominated era. There, logarithmic corrections or a different 
power law might have to be applied. Since the relevant scales enter the horizon roughly during the 
matter-radiation transition, this renders the amplitude of the corresponding fluctuations somewhat 
uncertain. 
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with particular solutions 

Dr = -—eAt^/^ ■ Dr = . ( 10 ) 

15 7 ^ ^ 

In the above equations we have introduced e = ATiCrf, the only free parameter in 

the model. We consider perturbations seeded by the texture held, and therefore 

it is incorrect to add a homogeneous growing mode to the above solutions. With 

these initial conditions, Eqs. (§), (|^) are easily integrated numerically, leading to 

the spectra for Dr{k,trec) and Dr{k,trec) [see, Fig. 1]. 



Figure 1: The dimensionless power spectra, k^\Dr\‘^ (solid line) and (dashed line) 

in units of {Ae)"^, are shown as functions of k. These are exactly the quantities which enter 
in the expression for the Cr’s. We set h = 0.5 , VLb = 0.05 and Zrec = 1100. 

Integrating Eq. (0), we obtain the Doppler contribution to the CMB anisotropies 
[see, Fig. 2]. For i < 1000, we hnd three peaks located at i = 365, i = 720 and 
i = 950. Silk damping, which we have not taken into account here, will substantially 
decrease the height of the second and even more that of the third peak. The inte¬ 
grated Sachs Wolfe effect, which also has been neglected, will shift the position of 
the hrst peak to somewhat larger scales, lowering ipeak by (5 - 10)% and increasing 
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Figure 2: The angular power spectrum for the Doppler contribution to the CMB 
anisotropies is shown in units of e^. We set cosmological parameters h = 0.5 , Qb = 0.05 
and Zrec = 1100 


its amplitude by less than 30%, as we argued above. 

Our second result regards the amplitude of the hrst Doppler peak, for which we 
hud 




^~365 


56 ^ . 
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As a consequence of the remark in footnote 4, the above numerical result has to be 
taken with a grain of salt. It is interesting to notice that the position of the hrst 
peak is displaced by Ai ~ 150 towards smaller angular scales than in inhationary 
models 0. This is due to the fact that our solution represents a combination of 
the growing and decaying modes, and only once the perturbation enters the horizon 
and the source term becomes negligible, the decaying mode starts to decay. This 
is manifest in the diherence in the growth of super-horizon perturbations, which is 
Dr oc in our case, and Dr oc for inhationary models, where on all scales only 
the growing mode is present. 

One may understand the height of the hrst peak from the following analytic 
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estimate: matching the sub-horizon with the super-horizon solutions of Eq. (^, in 
the matter dominated era, one hnds -Dc ~ —0. 474e(fc/27r)^/^t^. From Eq. (|^) we then 
obtain in this limit Dr ~ Aek~^^'^. Plugging this latter value into Eq. (i we get 
roughly i{i + l)Ci ~ for the height of the hrst peak. This agrees, within a 

factor 2 with the numerical result given in Eq. (0)- 

Let us now compare our value for the Doppler peak with the level of the SW 
plateau. Unfortunately, the numerical value for the SW amplitude is uncertain 
within a factor of about 2, which leads to a factor 4 uncertainty in the SW contri¬ 
bution to the power spectrum: Refs. 01 and Ref. [§] hnd respectively 


' + l)Ci 


SW 


2 e^ and £(£ + 1)Q 


SW 


8e^ 


( 12 ) 


According to Refs. [0 |^, the Doppler peak is a factor of ~ 3.4 times higher than 
the SW plateau, whereas it is only about 1.5 times higher if the result found in Q] 
is assumed. (We allow for about 30% of the SW amplitude to be added in phase to 
the Doppler amplitude of ~ 2.24e, according to Eq. O). Clearly, improved numer¬ 
ical simulations or analytical approximations are needed to resolve this discrepancy. 
However, it is apparent from Eqs. ([0]) and (^) that the Doppler contribution from 
textures is somewhat smaller than for generic inflationary models. 

We believe that our results are basically valid for all global defects. This depends 
crucially on the 1/v^ behavior of (00^ large scales (cf. Eq. (|^)), which is a 
generic feature of global defects: on super-horizon scales, (0')^(A;) represents white 
noise superimposed on the average given by = 0) oc \/V/D. Since there are 

N = V/t^ independent patches in a simulation volume V, the amplitude of 
is proportional to \/V oc (Notice that this argument does not apply 

for local cosmic strings.) 


Based on our analysis, we conclude that if the existence of Doppler peaks is in¬ 
deed conhrmed and if the hrst peak is positioned a.t £ < 300, then global topological 
defects are ruled out. On the other hand, if the hrst Doppler peak is positioned 





at £ ~ 350 and if its amplitude is lower than the one predicted for standard infla¬ 
tionary models, global topological defects are strongly favored if compared to the 
latter. (There are however non-generic, open, tilted inflationary models which might 
reproduce a similar signature in the CMB angular power spectrum). To our knowl¬ 
edge this is the hrst clear hngerprint within present observational capabilities, to 
distinguish among these two competing models of structure formation. 

As a future work, we aim to model with better accuracy the global scalar held 
0 during both the radiation and the matter dominated era, as well as to include 
the SW effect and the photon diffusion. This will allow us to better estimate the 
amplitude of the hrst Doppler peak and to investigate secondary peaks. 

As we were completing our work, a preprint]^ on the same issue, but follow¬ 
ing a diherent approach, came to our attention. The authors calculate the Doppler 
peaks from cosmic textures in the synchronous gauge. They include the Sachs Wolfe 
contribution into the analysis, but they need more of the uncertain modeling of the 
defect stress energy tensor. Even though we basically agree with the shape and 
position of their Doppler peaks, we obtain a somewhat smaller amplitude. 
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